Permanent polarization of small metallic particles 
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Electric charge density in a metallic particle fluctuates due to inhomogeneities of various kinds. 
While in the bulk of the particle the charge fluctuations are suppressed by Thomas-Fermi screening, 
the underscreened charges near the surface give rise to a permanent electric dipole moment. We 
show that the dipole moment increases linearly with the particle size, and fluctuates strongly from 
particle to particle. 
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Small metallic particles, even electrically neutral, gen- 
erate an electric field around them. The field arises be- 
cause the charge density within a particle is not uniform. 
Generally, the charge density includes both a dynamic 
(time-dependent) contribution, and a static, one. While 
the former is the origin of van der Waals interaction, 
the latter results in a formation of a finite static elec- 
tric dipole moment. The static dipoles may determine 
the interaction between particles at large distances. De- 
tailed understanding of the interaction between nanopar- 
ticles is important, e.g., for a self-assembly of functional 
nanostructures. 1 In particular, the dipolar interaction is 
responsible* for the tendency of nanoparticles to arrange 
in quasi-one-dimensional chains. 

The dipole moments of nanoparticles can be measured 
in Stern-Gerlach-type beam deflection experiments,^ or 
by studying the electric response of dilute solutions of 
nanoparticles, 5 Although none of these methods produces 
a reliable set of data for a sufficiently broad range of par- 
ticle sizes, measured values of the static dipole moment 
apparently increase with the particle size. 

While for very small particles (with number of atoms 
N < 10) the formation of the dipole moment can be 
studied quantitatively by ah initio methods^, the prob- 
lem becomes untractable for larger N . Indeed, the spatial 
arrangement of atoms may differ considerably from par- 
ticle to particle. The variations in shape and structural 
defects make the problem even more complex. 

One can view a metallic particle as consisting of a rigid 
"skeleton" formed by charged ions, and conduction elec- 
trons moving freely on this background. The skeleton's 
charge distribution includes a random component repre- 
senting structural defects, surface roughness, etc. Such 
structural disorder is present even if the bulk material 
with the same chemical composition is crystalline. It is 
therefore natural to analyze properties of nanoparticles 
statistically, averaging quantities of interest over the dis- 
order realizations. In this paper, we apply this strategy 
to the problem of evaluation of the static dipole moment. 
We compute the disorder average of p 2 and find 



<P 2 > 



e 2 S, 



(1) 



where S ~ L 2 is the surface area of the particle, and 



e is the charge of an electron. (Note that (p) vanishes 
on symmetry grounds.) We derive Eq. (JTJ) by consider- 
ing the disorder-induced mesoscopic density fluctuations. 
We show that the fluctuations originate on a short scale 
of the order of an interatomic distance. The density fluc- 
tuations are screened in the bulk of the particle. Near 
the particle's surface, however, the fluctuating charges 
and their screening "clouds" form dipoles. Adding these 
random dipoles, we arrive at Eq. |T]), which implies that 
the typical value of th e dip ole moment scales with the 
particle size L as p ~ yj (p 2 ) ~ eh. As discussed towards 
the end of the paper, the scaling p ~ eL is in apparent 
agreement with the existing experimental data. 

We start the analysis by noting that in a metallic parti- 
cle with a large number of electrons, all electronic states 
below the Fermi level ef contribute to the static elec- 
tron density n(r). Quantities of this type can usually be 
estimated semiclassically. Such a consideration yields a 
homogeneous density n(r) that compensates exactly the 
ion background charge, n(r) — n lon = const. However, 
the semiclassical description breaks down in the vicinity 
of defects or boundaries. Indeed, an impurity embedded 
in a Fermi gas perturbs the density around it. The excess 
density Sn(r) — n(r) — m oa is described by the Friedel 
oscillation which decays with the distance r to the impu- 
rity as (Xf/t) 3 , and oscillates with the spatial period of 
the Fermi wavelength \p. (Note that disorder does not 
destroy Friedel oscillations^) Accordingly, the disorder- 
induced density fluctuations are statistically correlated 
on the scale of Xf- 

In order to estimate the short-range contribution to 
the density fluctuations, we consider an infinite system 
with randomly located identical pointlike impurities. By 
the Friedel sum rule, the change of the total number 
of electrons due to a single impurity is given by Sq/tt, 
where Sq is s-wave scattering phase shift. To lowest order 
in | (So | <C 1, contributions to 5n(r) from different impuri- 
ties are additive. Neglecting the details of the charge dis- 
tribution on the scale of the order of Xf, one can relate 
the excess density to the local impurity concentration, 
5n(r) m 2(5q / TT)m(r) , with the factor of 2 accounting for 
the spin degeneracy. Averaging over the positions of the 



impurities then yields 

K(r - r') = (Sn(r)Sn{r')) = An, (S /tt) 2 S(r - r'), (2) 

where n, = (rii(r)) is the average impurity concentration. 
The phase shift So can be expressed via the impurity 
scattering cross-section Oi — X F S 2 /n. Relating the latter 
to the elastic mean free path I — (rijCTj) -1 , we rewrite 
Eq. |2]) as 



K(r - r') = K hulk S(r ~ r') , 
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-Kbulk = TT7 ■ (3) 
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Below we show that the estimate © coincides with 
the result of a rigorous calculation. Before doing so, 
we compare Eq. ([3]) with the results of Refs. [1| and 
Q. These authors considered only a contribution of 
electronic states near the Fermi level, treating them in 
the diffusion approximation. Due to a relatively small 
number of such states, the corresponding contribution 
to K (r — r') is smaller by a factor (Xp/l) 2 <C 1 than the 
estimate ©. For the same reason, a small factor (Xp/l) 2 
appears in the result of Ref. Q for the static dipole mo- 
ment. Thus, Eq. ([3]) indeed describes the dominant con- 
tribution to the static density fluctuations in a disordered 
metal. It should be emphasized that due to a large num- 
ber of electronic states contributing to density fluctua- 
tions, Eq. ([3j is not affected by quantum interference 
effects such as weak localization. 

We now calculate the static density-density correlation 
function microscopically. We consider first noninteract- 
ing electrons described by the single-particle Hamiltonian 



H 



2m 



U(r) 



(4) 



Since the density fluctuations are correlated on the scale 
Xf <C L, we ignore the boundaries for the time being and 
consider an infinite system. The random potential U (r) 
in Eq. (j3J) represents the structural disorder, which is 
correlated at distances of the order of the interatomic 
spacing. For simplicity, we assume this correlation length 
to be small compared to Xf, and write 

<[/(r)) =0, (U(r)U(r')) = -^-d(r-r'), (5) 

ZlTVT 

where v — mpp/^n 2 is the density of states, r is the 
mean free time, and pf is the Fermi momentum. (A more 
realistic form of Eq. ^ would only affect the numerical 
coefficient in the expression for ifbuik-) Introducing the 
causal Green's function G(i,r,r') = — i(TV> + (i, r )V'( r/ ))j 
and taking into account the spin degeneracy, we write 
the electron density as n(r) = 2iG(+0, r, r). In frequency 
domain, this relation takes the form 



n(r) 



de Im G(e, r, r) 
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FIG. 1: (a) The lowest order contribution to the static 
density-density correlation function. The solid and dashed 
lines represent electron Green's functions and disorder poten- 
tial, respectively, (b) Screening cloud in a metal. A point 
charge q in the bulk induces a symmetric screening cloud. 
A charge q' near the surface induces an asymmetric screening 
cloud with the screen charge —q centered at a different point. 
This offset results in a finite dipole moment. 



The lowest order in Xf/1 contribution to the product of 
two Green's function corresponds to the diagram-^ shown 
in Fig. [UJa) , which gives 

A ' (r - r '> -55? /s£/ dV '« r - -">«■ - r '> ■ 



Here F e (r) = ImGQ(e,r) with 



(7) 



Go(e,r) 



a ip e r sgn(e-e F ) 
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being Green's function in the absence of disorder. Phys- 
ically, Eq. ((?]) represents a product of two Friedel oscilla- 
tions induced at points r and r' by a single impurity, av- 
eraged over the impurity position r" . Analysis of Eq. ([7]) 
shows that K (r — r') changes on the scale of the order 
of Xf, which is the shortest length scale in our problem. 
At larger distances, K(r — r') can be approximated as 

K(r - r') a K hulk S(r - r') , K hulk = Jd 3 rK(r) . 

(8) 

Substituting Eq. ([7]) into Eq. ([8|) and using pf = 2-k/Xf 
and t = ml/pF, we obtain if bulk given by Eq. (|3|). 

Formally, Eq. © assumes a constant chemical poten- 
tial and is valid only for an infinite system. In a finite 
system, however, it is the number of electrons that is 
kept constant. A charge pushed away by an impurity is 
spread over the volume of the particle. This can be ac- 
counted for by introducing an appropriate correction to 
the chemical potential. For L ^ Xf, this amount*^ to 
the replacement 8(r — r') — * <5(r — r') — V^ 1 in Eq. ([H]); 
here V ~ L 3 is the particle volume. This yields 

K(t) w ifbuik [5{r) - 1/V) , [ d 3 rK(r) = 0, (9) 

Jv 

which is consistent with charge conservation. 

We turn now to the evaluation of the dipole moment. 
It is related to the density fluctuation 5n(r) as 



p = e J d 3 rv5n{r) 



(10) 
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The random vector p vanishes upon averaging over the 
disorder realizations, (p) = 0. Its dispersion (p 2 ) = (p 2 ), 
however, is finite and is expressed via the density-density 
correlation function, 

(p 2 ) = e 2 J d 3 r J d 3 r'(r ■ r')(Sn(r)Sn(r')) . (11) 

Using here Eq. © for (Sn(r)5n(r')) , one would find 

(p 2 ) cx e 2 K bulk L 5 , (12) 

i.e., p grows with particle size L as L 5 / 2 . However, Eq. (J9j) 
is valid only in the absence of Coulomb interaction be- 
tween electrons. The interaction suppresses charge fluc- 
tuations at distances larger than the screening length l s . 
Thus, Eq. (fT2")l is applicable only for L < l s . Since l s 
is rather short for most metals, this condition is very 
restrictive, and interaction must be taken into account. 

Treating disorder and interaction simultaneously is 
usually a very difficult task. It is greatly simplified when 
the screening length l s is large compared with the Fermi 
wavelength A^. In the limit l s 3> Xf, one can analyze 
the effects of the disorder and of the interaction indepen- 
dently of each other. At short distances, |r — r'\ <C l a , the 
screening is not effective, and the density fluctuations are 
still described by Eq. ([3]) . The fluctuations are "dressed" 
by the screening "clouds" at longer distances. The dress- 
ing can be described in Thomas-Fermi approximation, 
which relates the screening density 8n s (r) to the elec- 
trostatic potential <p(r) as 5n s (r) = —Kecp(r). Here k is 
the compressibility of the electron gas; disorder-induced 
fluctuations of k are negligible The potential ip(r) 
satisfies Poisson's equation 

W 2 ip = -47re<5n(r), (13) 

where the total density 5n(r) = Sno(r) + 5n s (r) includes 
both the bare density 5no(r) described by Eq. (jjJJ) and 
the screening contribution Sn s (r). In general, Thomas- 
Fermi approximation assumes a constant chemical poten- 
tial. To apply it to a finite size system, one has to in- 
troduce a correction to the chemical potential to ensure 
charge conservation. Such a modification is not needed 
in our case, since the total charge of the bare density 
fluctuations vanishes, J d 3 r 5no(r) = 0. Therefore, the 
total screening charge vanishes as well. 

Excluding (p(r) from Eq. HH), we write Sn(r) as 

5n(r) = 5n {r) - l~ 2 J dV Q(r, r') Sn (r') (14) 

where l s = (47re 2 K)~ 1 / 2 is the screening length, and the 
function Q(r, r') satisfies the equation 

[-V 2 + C 2 ] Q(rX) = <5(r - r'). (15) 

This equation has to be supplemented with the boundary 
condition which corresponds to vanishing of the electric 
field far away from the particle. That is, one has to solve 



Laplace equation V 2 Q = in the exterior of the particle 
and match the values of Q at the particle boundary. 

The screening described by Eq. (fT"5|) occurs on the scale 
l s . For a large particle, L ^> l s , this scale is small com- 
pared to the radius of curvature (~ L) of the particle's 
boundary, which allows us to treat the boundary as be- 
ing locally flat. We choose the coordinate axes so that 
x > corresponds to the interior of the particle. It 
is convenient to integrate Eq. (fT5|) along the boundary, 
which makes the electrostatic problem one dimensional. 
Eq. (fT5"l) preserves its form, except for V 2 being replaced 
by d 2 /dx 2 . Solution of the Laplace equation in the ex- 
terior (x < 0) reads Q(x,x') = const, which gives the 
boundary condition dQ/dx = at x — 0. Solution in 
the interior (x > 0) can then be found by the method of 
images, 

J Q(r,r')dy'dz' = l s (e^ x - x '^ + e -0H-»')/».) . ( 16 ) 
With the help of the relations 

Jd\'Q(v.v') = I 2 , (17) 

Jd 3 r'r' Q(r,r') - x l 2 s (x + l s e~ x/l ^ , (18) 

we see that the total charge of the distribution ([M]) is 
zero, while the dipole moment is finite and is given by 

5p = el s Jd 3 rxSn {r)e' x/h . (19) 

Thus, the screening transforms a charge q at distance 
x from the surface into a dipole with the dipole moment 
~ ql s exp(—x/l s ). Qualitatively, this means that a charge 
in the bulk (at x S> l s ) is screened by a spherically sym- 
metric cloud. However, if the charge is close to the sur- 
face (at x < l s ), the screening cloud is asymmetric, see 
Fig. Hlb) > gi ym g r i se to the dipole moment of the order 
of ql s . 

For a smooth boundary of an arbitrary shape the co- 
ordinates x and x' in Eqs. (fT!))) - (fl"5)) should be replaced 
by the distances from points r and r' to the boundary, 
and vector x by the inward unit normal to the boundary 
at the point closest to r. Using Eq. (|14p. we write the 
density-density correlation function in the form 

(Sn(r)Sn(r')) = J d 3 n d 3 r 2 Q(r, r^i^n - r 2 )Q(r 2 , r'). 

(20) 

Here Q(r, r') = Q(r, r') — l 2 S(r — r'), and K {v\ — r 2 ) given 
by Eq. @ describes fluctuations of Sno(r). Substituting 
Eq. fl2DJ| into Eq. {TDJ) and using the relations (T7]) and 
fT5|). we find 

(p 2 ) = e 2 l 2 J d 3 rd 3 r' K(r - r >) e -l*M+*(r')]/l* , ( 2 1) 

where x(r) is the distance from point r to the surface. 
Taking into account that L 3> l s , we obtain 

(p 2 ) = le 2 l 3 s K hulk S= 1 Lp I S , (22) 
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where S ~ L 2 is the surface area of the particle, 
and -Kbuik given by Eq. §5§ characterizes the bulk dis- 
order. 

The main contribution to the integral in Eq. (|2ip comes 
from a thin layer of thickness l s -C L near the surface. 
Therefore, the dipole moment is sensitive to the details 
of the surface structure, such as roughness, adsorbed im- 
purity atoms, etc. The effect of the surface disorder can 
be also treated by our method. Indeed, similarly to the 
bulk impurities, the surface defects induce short-range 
density fluctuations. With the screening taken into ac- 
count, the corresponding contribution to the dipole mo- 
ment has the form of Eq. I|22p with inbuilt replaced by an 
analogous quantity characterizing the surface disorder. 
Thus, accounting for the surface disorder does not affect 
the scaling (p 2 ) ~ e 2 S. 

The above consideration allows us to find not only the 
dispersion (p 2 ), but the entire distribution function of the 
dipole moment. Indeed, since the bare density fluctua- 
tions are correlated only at short distances, the dipole 
moment is a sum of a large number (~ Sl s /X F 3> 1) 
of independent random contributions. By the Central 
Limit Theorem, the distribution function of p is Gaus- 
sian, /(p) oc exp(— p 2 /2(p 2 )). Accordingly, the distribu- 
tion of the absolute value p = |p| of the dipole moment 
is given by 

/W«P 2 exp(-^y). (23) 

In principle, the statistics of p can be studied in Stern- 
Gerlach-type experiments, see, e.g., Ref. |3|. 

Formally, the analysis given in this paper is valid only 
when l s ,l 3> Xf- In order to estimate the upper bound 
to the permanent dipole moment, we consider the lim- 
iting case l s w I w Af, and find p ~ ey/S ~ eL. 
The dipole moment of this order of magnitude corre- 
sponds to the transferring of one electron across the par- 



ticle. This agrees with the existing experimental data. 
In Ref. 0, the dipole moment of niobium (Nb) clus- 
ters with N Ri 90 atoms was found to be of the order 
of 0.5 Debye (D) per atom. Assuming the cluster to be 
a sphere of radius a{?>/ AttN) 1 ^ , where a rj 3 A is the 
lattice constant for Nb, we estimate the dipole moment 
per atom to be p/N ri e^S/An/N w 0.5 D. Accord- 
ing to Ref. [5j , 3 — 5 nm - size semiconducting (such as 
CdSe) nanocrystals have permanent dipole moments in 
the range 50 — 100 D. Even though our theory was de- 
veloped for metallic clusters and cannot be directly ap- 
plied to semiconductors, the estimate based on Eq. (|2"2"]) 
gives p ~ 100 D. 

Although the static dipole moment p grows linearly 
with the particle size L, the electric field induced by the 
particle in its exterior, E ~ pjl? ~ e/L 2 , vanishes in the 
limit L — > oo. Also, the polarization (dipole moment per 
volume) p/V vanishes in the limit L — > oo, as expected 
for a normal metal. This behavior is to be contrasted 
with that of ferroelectrics where p/V — ► const in the 
thermodynamic limit. 

On the other hand, the main contribution to the 
disorder-induced density fluctuations that give rise to the 
finite dipole moment comes from a large number of the 
deep-lying electronic energy levels^ 2 - Therefore, we ex- 
pect the static dipole moment to be independent of tem- 
perature T as long as T is small compared with the Fermi 
energy. However, at T > e 2 /L, thermal fluctuations of 
the dipole moment exceed its static value. 

To conclude, in this paper we have identified the dom- 
inant contribution to the charge density fluctuations in a 
disordered metal and applied the result to the evaluation 
of the permanent dipole moment of a metallic particle. 

We thank P. Allen, W. de Heer, and U. Landman for 
valuable discussions. This work is supported by the DOE 
(contract No. DEAC 02-98 CH 10886), and by the NSF 
(grant No. DMR-0604107). 



1 J.V. Barth, G. Costantini, and K. Kern, Nature 437, 671 
(2005); G.M. Whitesides and B. Grzybowski, Science 295, 
2418 (2002); F.X. Redl, K.-S. Cho, C.B. Murray, and S. 
O'Brien, Nature 423, 968 (2003). 

2 K.-S. Cho, D.V. Talapin, W. Gaschler, and C.B. Murray, 
J. Am. Chem. Soc. 127, 7140 (2005). 

3 R. Moro, X. Xu, S. Yin, and W.A. de Heer, Science 300, 
1265 (2003); R. Moro, S. Yin, X. Xu, and W.A. de Heer, 
Phys. Rev. Lett. 93, 086803 (2004). 

4 R. Schafer, S. Schlecht, J. Woenckhaus, and J. A. Becker, 
Phys. Rev. Lett. 76, 471 (1996). 

5 M. Shim and P. Guyot-Sionnest, J. Chem. Phys. Ill, 6955 
(1999); L. Li and A. P. Alivisatos, Phys. Rev. Lett. 90, 
097402 (2003). 

6 K.E. Andersen, V. Kumar, Y. Kawazoe, and W.E. Pick- 
ett, Phys. Rev. Lett. 93, 246105 (2004); Phys. Rev. B 73, 
125418 (2006). 

7 AYu. Zyuzin and B.Z. Spivak, JETP Lett. 43, 234 (1986) 
[Pis'ma v Zh.E.T.F. 43, 185 (1986)]; L.N. Bulaevskii and 



S.V. Panyukov, JETP Lett. 43, 240 (1986) [Pis'ma v 
Zh.E.T.F. 43, 190 (1986)]. 

8 R. Berkovits and B.L. Altshuler, Europhys. Lett. 19, 115 
(1992); Phys. Rev. B 46, 12526 (1992). 

9 S. Ulreich and W. Zwerger, Zeits. Phys. B 97, 151 (1995). 

A. A. Abrikosov, LP. Gor'kov, and I.Ye. Dzyaloshinskii, 
Methods of Quantum Field Theory in Statistical Physics 
(Dover, 1975). 

1 Ya.M. Blanter and A.D. Mirlin, Phys. Rev. B 57, 4566 
(1998). 

2 Unlike the permanent dipole moment, the dipole moment 
induced by an external field is due to states close to the 
Fermi level. This makes it possible to analyze such linear 
response using the diffusion approximation, see Refs @i, 

3 K.B. Efetov, Super-symmetry in Disorder and Chaos (Cam- 
bridge University Press, 1997), Ch. 7. 



